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摘要 ”中国 散 裂 中 子 源 (CSNS) 可 提供 0.3 eV 至 300 MeV 的 白光 中 子 束 流 ， 总 束 流 强 度 可 达 107n/s/cm2， 为 
中 子 俘获 反应 截面 的 测量 提供 了 一 个 优秀 的 实验 平台 。CSNS 在 正常 的 运行 模式 下 会 由 加 速 器 产生 两 个 间隔 

局 先后 打靶 ， 因 此 产生 的 中 子 束 流 也 由 间隔 为 410 ns 的 两 个 束 团 混合 而 成 。 为 避免 两 个 束 团 
的 效应 相互 干扰 ， 影 响 中 子 爷 获 截面 的 能 量 精度 ， 需 要 对 实验 数据 进行 解析 和 重 构 ， 还 原单 个 束 团 的 效应 。 
见 有 的 解析 方法 可 以 得 到 非常 精细 的 解 谱 结果 ， 但 相对 复杂 ， 具 有 一 定 使 用 门槛 。 本 工作 提出 了 一 种 简化 的 
双 束 团 解 谱 方法 ， 在 保证 中 子 能 量 精度 的 情况 下 适用 于 中 子 能 量 低 于 1.2 MeV 的 数据 ， 为 同类 型 的 实验 工作 
提供 一 种 新 的 数据 处 理 思路 。 
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Abstract [Background]: The China Spallation Neutron Source (CSNS) provides a white neutron beam with an 
energy range from 0.5 eV to 300 MeV and a total beam intensity of up to 10’n/s/cm’, serving as an excellent 
experimental platform for the measurement of neutron capture reaction cross sections. During normal operation, the 
CSNS generates two proton bunches separated by 410 ns, consecutively striking the target, resulting in a mixed 
neutron beam composed of two bunches with a 410 ns interval. To avoid interference between the effects of the two 
bunches and maintain the energy precision of neutron capture cross sections, experimental data need to be analyzed 
and reconstructed to restore the effects of individual bunches. [Purpose]: The existing parsing method can yield very 
refined unfolding results, but it is relatively complex and has a certain usage threshold. Therefore, a more convenient 
data processing method needs to be found. [Methods]: This work utilized mathematical operations to analyze and 
reconstruct the data, with 410 ns as the unit time, and processed the data with a channel width of 4100 ns. 
Additionally, a comparison was made of the impacts of this method and existing methods on the accuracy of neutron 
incident energy. [Results]: This work proposes a simplified data processing method that achieves the same energy 
resolution as existing methods in the low-to-medium energy range, providing a new data processing approach for 


similar experimental work. [Conclusions]: The simplified data processing method presented in this study effectively 
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addresses the issue of excessive computational costs in analyzing low to medium energy neutron data from the CSNS. 
It offers a practical solution for experimental work requiring accurate analysis of neutron capture reactions in this 
energy range. 


Key words CSNS,Double-bunch unfolding method, Neutron capture reactions 


中 国 散 裂 中 子 源 (CSNS) 位 于 广东 省 东莞 市 ， 是 我 国 第 一 个 白光 中 子 源 实验 平台 ， 为 中 
国力 至 全 球 的 科研 工作 者 提供 了 一 个 强大 的 研究 平台 02。CSNS 的 中 子 束 流 是 利用 加 速 器 加 
速 并 经 过 磁场 偏转 后 的 高 能 质子 束 流 友 击 散 裂 靶 产生， 其 加 速 器 频率 为 25 Hz， 质 子 能 量 约 
为 1.6 GeV。CSNS 装置 概况 如 图 1 所 示 ， 与 质子 束 流 打 部 方向 呈 180? 的 角度 引出 一 条 中 子 
束 流 ， 称 之 为 反 角 白光 中 子 束 流 线 (back n) ， 包含 了 从 0.3 eV 至 300 MeV 范围 内 能 量 连续 
的 中 子 ,最 大 中 子 通 量 可 达 107n/cm2ks。 距 离 打 靶 点 55 米 和 76 米 处 分 别 设置 有 ES1#1 和 ES#2 
两 个 实验 终端 ， 其 中 ES1#1 具有 更 高 的 中 子 束 流 强度 ， 而 ES#2 具有 更 高 的 中 子 能 量 分 辨 。 
根据 需求 ， 用 户 可 以 选择 合适 的 终端 开展 实验 B61。 自 2018 年 打靶 出 束 以 来 ， 为 基础 物理 、 
材料 科学 等 多 个 学 科 的 实验 提供 了 束 流 支 持 5-9。 


1.6 GeV protons 


1 CSNS 装置 示意 图 ， 引 自 文献 [9]。 
Fig.l Schematic Diagram of the CSNS Facility, cited from Refl9]. 

在 常规 运行 时 ， 加 速 器 采用 的 是 双 束 团 模 式 ， 即 两 个 相同 的 质子 束 团 以 410 ns 为 间隔 裴 
击 散 裂 靶 。 因 此 在 实验 中 测量 到 的 是 两 个 间隔 为 410 ns 的 反应 效应 的 县 加 。 实 验 中 的 中 子 入 
射 能 是 根据 飞行 时 间 计 算 的 ， 时 间 的 起 点 为 第 一 个 质子 束 团 的 打靶 时 间 ， 第 二 个 中 子 束 团 中 
相同 能 量 的 中 子 对 应 的 飞行 时 间 则 对 应 增加 了 410ns。 在 测量 中 子 俘获 反应 截面 时 ， 两 个 束 
团 的 登 加 会 导致 计算 出 的 中 子 的 入 射 能 精度 下 降 ， 对 于 中 高 能 的 中 子 影响 尤为 显 普 。 因 此 ， 
为 了 得 到 准确 的 中 子 俘获 反应 激发 函数 ， 双 束 团 的 数据 必须 进行 解析 重 构 ， 还 原 出 单个 个 中 
子 束 团 真 实 的 效应 。 

中 科 院 高 能 物理 研究 所 的 易 蛤 等 人 在 2020 年 针对 CSNS 双 束 团 数据 开发 了 一 套 解析 程序 
0 ， 可 以 得 到 精细 的 中 子 反 应 截面 ， 用 户 可 以 在 此 网 址 下 载 编译 
(http://code.ihep.ac.cn/yih/csns-back-n-doublebunchunfolder〉 。 使 用 时 需要 先 将 测量 得 到 的 数 
据 转 化 成 直方 图 样式 的 root 文件 ， 再 将 该 程序 租 入 到 数据 处 理 程序 中 ， 经 过 30 次 至 上 百 次 
的 反复 迭代 ， 最 终 得 到 飞行 时 间 单 的 一 维 谱 或 者 飞行 时 间 与 能 量 的 二 维 谱 。 为 保证 解 谱 后 的 
数据 能 还 原 真实 的 反应 产 额 ， 推 荐 直方 图 的 时 间 道 宽 不 超过 延 时 的 十 分 之 一 ， 即 41 ns。 传 统 
的 中 子 俘获 截面 的 测量 使 用 CeDs 液 体 闪 烁 探测 器 等 具有 较 高 伽 马 射 线 灵 敏 度 的 探测 右 , 测量 
过 程 中 获取 反应 中 产生 的 所 有 人 徊 马 射 线 ， 数 据 处 理 时 通过 权重 因子 计算 出 中 子 俘获 反应 的 比 
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例 ， 进 而 得 到 目标 反应 的 截面 。 因 此 探测 器 获取 的 统计 量 通 常 足够 大 ， 即 使 将 时 间 谱 的 道 宽 
取 非 常 小 的 数值 ， 也 可 保证 较 小 的 统计 误差 。 这 样 的 处 理 方法 可 以 得 到 非常 精细 的 激发 函数 
曲线 《〈 中 子 俘获 反应 截面 随 入 射 能 变化 的 曲线 ) ， 但 也 会 带 来 巨大 的 计算 量 ， 因此 程序 的 作 
者 也 建议 数据 解析 的 时 间 跨 度 不 超过 105 ns， 否 则 可 能 会 导致 电脑 内 存 错误 ,程序 无 法 运行 。 
此 外 , 该 程序 的 基于 Linux 系统 开发 , 且 使 用 方法 相对 复杂 , 具有 一 定 的 门槛 , 对 不 熟悉 Linux 
系统 以 及 root 程序 代码 操作 的 用 户 有 一 定 难度 。 

对 于 两 种 情形 : 1) 实验 数据 的 统计 量 不 足 ， 使 用 较 小 的 时 间 道 宽 意味 着 每 道 的 统计 量 很 
小 ， 会 导致 统计 误差 过 大 ; 2) 测量 低能 中 子 的 数据 ， 即 飞行 时 间 在 105 ns 至 106ns 区 间 ， 宽 
道 宽 对 能 量 误差 的 影响 变 得 可 以 忽略 ， 而 继 纪 卖 使 用 较 小 的 时 间 道 宽 意味 着 极 大 的 运算 量 ， 不 
具有 必要 性 。 因 此 ， 上 述 两 种 情形 更 适合 采用 较 大 的 时 间 道 宽 ， 此 时 选择 简化 的 解 谱 方法 可 
以 在 符合 实验 精度 的 前 提 下 ， 实 现 数 据 解析 的 目的 。 本 工作 以 4100 ns 作为 时 间 窗 口 为 例 ， 
提出 了 一 种 双 束 团 数据 的 简化 解析 方法 ， 为 同类 型 的 中 子 俘获 反应 截面 测量 提供 一 种 新 的 数 
据 处 理 方案 。 


1 双 束 团 数据 解析 原理 


Back_n 的 中 子 束 流 可 以 看 做 是 两 个 能 量 分 布 相 同 ,但 是 前 后 相差 410 ns 的 中 子 束 团 登 加 ， 
分 别称 为 束 团 A 和 束 团 B。 探测 器 的 时 间 以 束 团 A 产生 的 时 刻 为 起 点 。 因 此 ， 在 探测 器 时 间 
0~ 410 ns 内 , 探测 器 获取 的 数据 全 部 来 自 于 束 团 A 在 这 个 时 间 段 的 效应 ; 在 探测 器 时 间 410 
~ 820 ns 内 ， 则 包括 了 束 团 A 在 第 二 个 时 间 段 的 数据 以 及 束 团 B 在 第 一 个 时 间 段 内 的 数据 。 
依次 类 推 , 在 之 后 的 每 个 410 ns 时 间 段 内 , 均 包 括 了 束 团 A 在 该 时 间 段 内 的 数据 以 及 束 团 B 
在 上 一 个 时 间 段 内 的 数据 。 以 410 ns 为 单位 时 间 ， 假 设 单个 束 团 在 第 n 个 单位 时 间 段 内 产生 
的 数据 为 fm， 探测 器 在 对 应 时 间 段 内 获取 的 数据 为 Am)。 那 么 有 : 

A(1)= {£1) 

A(2)=f(1)+ f(2) 


A(n+1) = f(tn)+ fln+1) (1) 

由 于 时 间 窗 口 为 4100 ns， 每 组 数据 包括 了 10 个 时 间 段 。 令 其 中 任意 一 组 数据 的 时 间 起 

点 为 第 x 个 时 间 段 的 起 上 点。 单个 束 团 在 该 时 间 窗 口内 的 数据 之 和 S=f(x)+f(x+1) 

Hf(x+2)+...+f(x+9)。 由 式 (1) 可 知 ，flx)+f(x+1)=A(x+1)，fl(x+2)+f(x+3)=A(x+3)，... 

f(x+8)+f(x+9)=A(x+9)。 因 此 单个 束 团 的 数据 与 探测 器 获取 的 数据 之 间 的 关系 式 可 以 改写 为 : 
S=A(x+1) + A(x+3) + A(x+5) + A(x+7) + A(x+9) (2) 

由 式 (2) 可 知 , 间隔 取出 数据 再 重新 组 合 就 可 以 得 到 该 组 数据 中 单个 束 团 的 数据 。 此 外 ， 

由 于 探测 器 对 于 不 同 能 量 的 伽 马 射线 具有 不 同 的 探测 效率 ,还 需要 取出 每 组 数据 的 能 量 单 谱 。 
本 工作 中 将 数据 解析 过 程 分 为 四 个 步骤 : 1) 将 存储 在 root 文件 中 的 时 间 与 能 量 二 维 谱 导 

出 为 文本 文档 2) 以 410 ns 为 时 间 间 隔 将 文档 拆 分 成 多 个 子 文档 ，3) 对 每 组 数据 间隔 取出 

子 文档 ， 合 并 成 一 个 新 的 文档 ，4) 获取 每 组 数据 的 能 量 单 谱 。 
此 外 ， 由 式 (1) 和 式 (2) 可 知 ， 本 方法 最 小 的 时 间 道 宽 为 820 ns。 


2” 双 束 团 数据 解析 过 程 


2.1 ”导出 root 数据 为 文本 文档 


日 于 不 同 数据 获取 系统 生成 的 root Branch 结构 , 此 处 以 本 团队 的 数据 结构 
为 例 。 其 中 ， 能 量 和 时 间 的 Branch 名 称 分 别 为 E 和 DeltaT， 数 据 类 型 分 别 是 Double t 和 
Oona 分 别 定义 两 个 对 应 数据 类 ee 


EH 


XXXXXX—3 


核 技 术 20XX，XX: XXXXXX 


间 数 据 。 完 整 的 数据 处 理 代码 见 附录 1。 代 码 的 思路 是 逐 事件 读 取 root 文件 数据 ， 并 将 每 个 
事件 的 能 量 与 时 间 信 息 存储 到 文本 文档 中 ， 用 户 在 使 用 时 根据 具体 情况 修改 对 应 内 容 。 以 下 
为 关键 步骤 的 程序 语句 : 

1) 实现 变量 与 数据 的 对 应 : 

tl->SetBranchAddress("E", &E1); Vtl 是 Tree 的 指针 。 

tl->SetBranchAddress("DeltaT", &DT1); 

2) 创建 文本 文档 input.dat， 用 以 存储 导出 的 数据 : 

ofstream outfile("input.dat"); 

3) 获取 root 文件 的 总 事件 数 ， 并 逐 事件 运行 : 

Long64 tntotl = tl->GetEntriesFast(); // 获 取 总 事件 数 

for(i= 0; 1< ntotl; 1++) 


{ 


tl1->GetEntry(i); 


outfile << DT1 <<" "<< El1<<endl; 


} 


2.2 ” 拆 分 数据 文档 

本 部 分 完整 代码 见 附录 2。 程 序 运行 时 需要 读 取 intervals.txt 文件 ， 其 中 存储 了 时 间 间 隔 
的 数值 ， 每 两 个 数值 为 一 组 ， 依 次 为 开始 和 结束 时 间 ， 单 位 为 ns。 拆 分 好 的 子 文件 依次 以 
001.dat、002.dat... 进 行 命名 ， 文 件数 量 由 时 间 的 组 数 决定 。 为 确保 遵循 数据 处 理 原 理 ， 开 始 
时 间 必 须 为 410 ns 的 整数 倍 ， 并 且 每 组 的 时 间 间 隔 均 为 410 ns。 


2.3 ”重组 数据 文档 


i 程序 运行 时 需要 读 取 datfile.txt 文件 ， 其 中 存储 了 所 有 子 文件 
的 文件 名 ， 包 含 后 级 。 重 组 后 的 数据 文件 依次 以 add1.dat、add2.dat... 进 行 命名 ， 文 件数 量 由 
0 程序 中 设 定 的 数据 组 上 限 为 44， 如 果 数 据 组 的 数量 超过 4$， 需 要 修改 程 
序 34 行 中 group 的 上 限 。 

2.4 ”获取 能 量 单 谱 

在 2.3 节 新 生成 的 数据 包括 了 时 间 及 能 量 二 维 信息 ， 其 中 第 一 列 为 时 间 ， 第 二 列 为 能 量 。 
获取 能 量 单 谱 的 方式 为 读 取 整个 文件 中 的 所 有 能 量 数值 ， 并 以 用 户 规定 的 道 宽 统计 每 个 道中 
能 量 数值 出 现 的 次 数 。 完 整 的 程序 见 附录 4， 需 要 提取 能 量 单 谱 的 文件 名 存储 于 filename.txt 
文件 中 。 
改 为 .txt 的 新 文件 。 新 文件 中 有 两 列 数 据 ， 第 一 列 为 能 量 ， 第 二 列 为 对 应 的 统计 。 至 此 ， 用 
户 可 获 得 单个 束 团 在 多 个 时 间 窗口 内 的 能 量 单 谱 ， 


2.5 ”开发 环境 


操作 系统 : Windows 10 专业 版 
Root 版 本 : root V6.30/04 
C++ 代 码 编 译 器 : Visual Studio 2021 
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3 结果 与 讨论 


中 子 的 入 射 能 与 飞行 时 间 的 平 发 旺 反比 ， 在 时 间 道 宽 恒 定 的 情况 下 ， 随 着 中 子 入 射 能 的 
增 大 ， 中 子 能 量 的 相对 误差 也 呈 平 方 关系 增 大 。 在 ES#2 终端 测量 中 子 俘获 反应 截面 时 ， 对 
于 10 keV 以 上 的 能 量 区 间 采 用 41 ns 为 时 间 道 宽 ， 可 以 将 中 子 能 量 的 相对 误差 控制 在 15% 以 
内 。 本 方法 主要 针对 10 keV 以 下 的 能 量 区 间 ， 采 用 4100 ns 为 时 间 道 宽 也 可 以 实现 相同 能 量 
误差 的 控制 。 图 2 分 别 展示 了 以 4100 ns 和 41 ns 为 道 宽 时 的 中 子 能 量 误差 曲线 。 其 中 实 线 表 
示 4100 ns 道 宽 ， 点 虚线 表示 41 ns 道 宽 。 可 以 看 出 ， 两 种 道 宽 在 各 自 对 应 的 能 区 内 具有 近乎 
相同 的 相对 误差 曲线 。 


一 一 This work 
= 训 1 


Error (%) 


102 10% 10 10° WW! 10 WY 10% 1 1 
E neutron (keV) 


图 2 中 子 能 量 相对 误差 曲线 。 实 线 为 4100 ns 道 宽 的 结果 ， 点 虚线 表示 41 ns 道 宽 的 结果 。 
Fig.2 Neutron energy relative error curve. 

为 了 更 好 地 展示 中 子 能 量 误差 的 影响 ,本 团队 在 back_n 束 流 线 上 对 ”Zr 的 中 子 俘 获 截 面 
进行 了 测量 ， 时 间 范 围 从 11480 ns 至 191880 ns， 时 间 道 宽 为 4100 ns， 共 44 个 数据 点 ， 对 应 
中 子 的 能 量 范围 在 0.82 ~ 229 keV。 利 用 本 方法 获得 的 中 子 截面 如 图 3 所 示 ， 从 图 中 的 横 轴 误 
差 棒 长 度 可 以 看 出 ， 当 中 子 入 射 能 超过 10 keV 时 ， 中 子 能 量 的 误差 十 分 明显 ， 相 对 误差 最 高 
可 达 90% 左 右 ， 中 子 俘获 截面 的 变化 趋势 也 相对 平缓 ， 难 以 判断 是 截面 本 身 趋 于 稳定 还 是 因 
为 截面 被 平均 后 变 得 平缓， 核反应 截面 的 结构 信息 被 掩盖 了 。 当 入 财 能 逐渐 减 小 时 ， 中 子 的 
能 量 误 差 已 控制 到 很 小 的 数值 ， 测 量 到 的 中 子 俘获 截面 震荡 变 得 十 分 剧烈 ， 很 好 地 反映 了 该 
能 量 区 间 内 的 俘获 截面 共振 特性 。 
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3 利用 本 方法 得 到 的 ”Zr 中 子 俘获 截面 。 带 XX 和 YY 四 点 为 截面 数据 。 曲 线 为 简单 的 连 线 ， 方 便 描 述 截面 的 变 
Fig.3 The Average cross sections of the neutron capture reaction of ?1ZT obtained using the proposed method in this work. 
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由 此 可 知 ， 当 中 子 入 射 能 超过 10 keV 且 统 计 足 够 大 时 ， 使 用 精细 的 解 谱 程序 可 以 获得 更 


加 精确 的 俘获 截面 ， 能 准确 体现 反应 截面 的 结构 信息 ， 而 较 宽 的 时 间 道 宽 会 丢失 这 些 结构 信 
息 。 当 中 子 入 射 能 低 于 10 keV 时 ， 使 用 本 方法 同样 也 可 以 获得 十 分 精确 的 激发 函数 曲线 。 


值得 注意 的 是 ， 本 方法 可 支持 的 最 小 时 间 道 宽 为 820 ns， 在 保证 中 子 能 量 分 辩 在 15% 以 


内 ， 最 高 可 适用 于 中 子 入 射 能 为 1.2 MeV 的 数据 。 


4 结语 


中 国 散 裂 中 子 源 back_n 束 流 线 为 中 子 俘获 反应 的 截面 测量 提供 了 一 个 绝 佳 的 实验 平台 。 


但 是 在 常规 模式 下 中 子 束 流 是 由 两 个 间隔 为 410 ns 的 束 团 混合 而 成 的 。 为 了 得 到 准确 的 中 子 
入 射 能 ， 需 要 对 双 束 团 数据 进行 解 谱 重 构 。 虽 然 已 经 存在 一 套 精细 的 解 谱 方法 ， 但 是 使 用 相 
对 复杂 ， 用 户 使 用 具有 一 定 的 门槛 。 而 且 对 于 需要 使 用 较 大 时 间 道 宽 的 数据 ， 精 细 解 谱 方法 
的 必要 性 不 足 。 本 工作 根据 基础 的 数学 运算 ， 提 出 了 一 种 十 分 简化 的 解 谱 方法 ， 以 较 大 的 时 
间 道 宽 对 数据 进行 拆 分 再 重 构 。 所 有 程序 均 可 直接 在 Windows 系统 下 编译 运行 ， 简 单 易 用 。 


当 使 用 最 小 时 间 道 宽 820 ns 时 ， 此 方法 在 中 子 入 射 能 低 于 1.2 MeV 能 区 内 得 到 的 中 子 能 


量 分 辨 可 达到 精细 解 谱 法 的 水 平 ， 并 且 计 算 量 得 到 了 极 大 的 简化 。 忆 一 方面 ， 相 对 精细 解 谱 
方法 ， 增 大 了 时 间 道 帘 ， 也 可 以 提高 每 道 的 数据 统计 量 ， 减 小 统计 误差 ， 对 于 统计 不 足 的 实 
验 数据 具有 显著 的 优势 。 一 言 蔽 之 ， 对 统计 量 足 够 大 且 中 子 入 射 能 较 高 的 数据 ， 采 用 现 有 的 


解析 方法 更 具有 优势 。 而 对 于 统计 量 较 小 以 及 中 子 入 射 能 低 于 中 低能 区 的 数据 ， 采 用 本 方法 


具有 计算 量 小 ， 操 作 简 单 等 优势 。 
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附录 1 


#include <fstream> 


— 


2  #include <string> 

3 ”void ana_extract(const char *fin1)/ 注 意 ; 函数 名 必须 与 .C 的 文件 名 相同 ， 和 否则 报错 

4 { 

5 Double t E1;// 能 量变 量 

6 Long64 ti; // 事 件数 变量 

7 ULong64_t DT1; /飞行 时 间 ， 单 位 为 ns 

8 TFile* fRunl = new TFile(fin1); 

9 TTree* tl = (TTree*)fRun1->Get("tr");，// 创 建 Tree 指针 变量 t1， 获 取 root 文件 中 的 TTree 
10 tl->SetBranchAddress("E", &E!1); 

11 tl->SetBranchAddress("DeltaT", &DT1); 

12 Long64 tntotl = tl->GetEntriesFast0; /获取 总 事件 数 

13 ofstream outfile("input.dat); /创建 数据 存储 文档 

14 outfile << "DT" <<" " << "E" <<endl// 在 文档 的 第 一 行 输出 DT 和 忆 
15 for(i = 0; i< ntotl; it+) // 逐 事件 处 理 数据 

16 { 

17 tl->GetEntry(i); 

18 if(i%(ntot1/10)==0) std::cout<<"processing:"<<i*10/(ntot1/10)<<"%"<<std::endl; /显示 数据 处 理 进 度 ， 可 删除 
19 outfile << DT1 <<" "<< El <<endl; /将 时 间 和 能 量 数据 写 入 文档 

20 } 

21 fRun1->Close(;// 关 闭 root 文件 ， 释 放 内 存 

22 outfile0.close(); 

23 } 

附录 2 


1 #include <iostream> 


2 #include <fstream> 
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#include <sstream> 
#include <vector> 
#include <string> 
#include <iomanip> 

// 创 建 数据 的 结构 变量 
struct Data { 


double DT:; 
double E; 


int main() { 


std::ifstream intervalsFile("intervals.txt"); /打开 存储 时 间 间 隔 的 文件 


if (lintervalsFile.is_open()) { 
std::cerr << "Failed to open intervals.txt" << std::endl; 
return 1; 

} 

std::vector<std::pair<double, double> > intervals; 

double start, end; 

// 读 取 时 间 的 起 点 和 终点 

while (intervalsFile >> start >> end) { 
intervals.push back(std::make pair(start,end)); 

} 

intervalsFile.close(); 

std::ifstream dataFile("input.dat"); 

if (ldataFile.is_open()) { 
std::cerr << "Failed to open data.dat" << std::end!l; 
return 1; 

} 

std::string title; 

std::getline(dataFile, title); 

std::vector<Data> dataVec; 

double dt, e; 

while (dataFile >> dt >>e) { 
Data data; 
data.DT = dt; 
data.E = e; 


dataVec.push_ back(data); 
} 


dataFile.close(); 


int fileCount = 0; 
for (const auto& interval : intervals) { 
std::vector<Data> filteredData; 


for (const auto& data : dataVec) { 
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47 if (data.DT >= interval.first && data.DT <= interval.second) { 

48 filteredData.push_ back(data); 

49 和 

50 } 

51 std::stringstream filename; /创建 存储 每 个 时 间 区 间 内 数据 的 子 文档 

92 filename << std::setfill('0') << std::setw(3) << ++fileCount << ".dat"; 

53 std::ofstream outFile(filename.str()); 

54 if (loutFile.is_open()) { 

55 std::cerr << "Failed to create output file: " << filename.str() << std::endl; 
56 return 1; 

57 } 

58 for (const auto& data : filteredData) { 

59 outFile << data.DT <<"" << (int)(data.E + 0.5) << std::endl; /四舍五入 将 能 量 转 换 为 整数 
60 } 

61 outFile.close(); 

62 } 

63 std::cout << "Data filtering completed successfully." << std::endl; 

64 return 0; 

65 } 


附录 3 


1 #include <iostream> 

2  #include <fstream> 

3 ”#include <vector> 

4  #include <string> 

5 ” #include <sstream> 

6 /W 读 取 文件 名 列表 

《7 std::vector<std::string> ReadFileNames() 
3 

9 std::ifstream inputFile("datfile.txt"); 
10 std::vector<std::string> filenames; 

11 std::string filename; 

12 while (std::getline(inputFile, filename)) 
13 { 

14 filenames.push_ back(filename); 
15 } 

16 return filenames; 

1 于 


18 / 将 源 文件 的 内 容 追 加 到 目标 文件 
19 void AppendFileContent(const std::string& srcFile, const std::string& destFile) 


20 { 

21 std::ifstream src(srcFile, std::ios::binary); 

22 std::ofstream dest(destFile, std::i0s::binary | std::ios::app); 
23 if (!src || !dest) 


XXXXXX-9 


核 技 术 20XX，XX: XXXXXX 


24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
45 
46 
47 
48 
49 
50 
51 
52 
53 


{ 
std::cerr << "Error opening files!" << std::endl; 
return; 
} 
dest << src.rdbuf\); 
} 
int main() 
{ 
std::vector<std::string> filenames = ReadFileNames(); 
for (int group = 0; group < 45; ++group) /根据 数据 组 的 数量 修改 group 的 上 限 
{ 
std::string destFilename = "add"; 
destFilename.append(std::to_string(group + 1)); 
for (int 1= 3 - destFilename.length(); 1> 0; --1) 
{ 
destFilename.insert(destFilename.begin(), '0"); 
} 
destFilename += ".dat"; /添加 文件 后 缀 名 
std::ofstream dest(destFilename, std::ios::binary); 
for (inti= group * 10; i< std::min((group + 1) * 10, static_cast<int>(filenames.size())); +HD) 
{ 
让 (i %2 1!=0)/W 间隔 取 子 文件 
{ 
AppendFileContent(filenames[i], destFilename); 
} 
} 
} 
return 0; 
} 


附录 4 


MD oo I oO 信人 DD 一 


hk fk 
je 


#include <iostream> 
#include <fstream> 
#include <sstream> 
#include <map> 
#include <string> 
#include <vector> 
int main() { 
/ 输入 文件 流 
std::ifstream filenameFile("filename.txt"); 
if (!filenameFile.is_open()) { 
std::cerr << "Failed to open filename.txt" << std::endl; 


XXXXXX-10 


作者 甲 等 ， 中 文 题目 


12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
45 
46 
47 
48 
49 
50 


return 1; 


} 
/ 输出 文件 流 


std::vector<std::string> filenames; 


std::string filename; 
while (std::getline(filenameFile, filename)) { 
filenames.push back(filename); 


} 


filenameFile.close(); 
/ 统计 每 个 间隔 的 计数 


for (const auto& filename : filenames) { 


std::ifstream file(filename); 


if (!file.is open()) { 


std::cerr << "Failed to open file: " << filename << std::endl; 


continue; 
} 
std::map<int, int> count Map; 
int dt, e; 


int bin = 20; // 能 量 道 宽 


while (file >> dt >> e) { 
int interval = (e /bin) * bin; // 取 bin 后 
countMap[interval]++; 
} 
file.close(); 
// 创 建新 文件 


std::string outputFilename = filename.substr(0, filename.find last_of(".")) + ".txt"; 


std::ofstream outputFile(outputFilename); 

if (loutputFile.is_open()) { 
std::cerr << "Failed to create output file: 
continue; 


} 


for (const auto& pair : countMap) { 


outputFile << pair.first << " " << pair.second << std::endl; 


} 


outputFile.close(); 


} 


std::cout << "Data counting completed successfully." << std::endl; 


return 0; 


" << outputFilename << std::endl; 
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